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Spin-polarised radio-frequency currents, whose frequency is equal to that of the gyrotropic mode, will cause 
an excitation of the core of a magnetic vortex confined in a magnetic tunnel junction. When the excitation radius 
of the vortex core is greater than that of the junction radius, vortex core expulsion is observed, leading to a large 
change in resistance, as the layer enters a predominantly uniform magnetisation state. Unlike the conventional 
spin-torque diode effect, this highly tunable resonant effect will generate a voltage which does not decrease 
as a function of rf power, and has the potential to form the basis of a new generation of tunable nanoscale 
radio-frequency detectors. 
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The physics of spin-transfer torque phenomena, and most notably spin-transfer induced high frequency dynamics, has recently 
become a field of much experimental interest, where motion and oscillations of the magnetisation in nanostructures can be 
induced by the application of a spin-polarised current. These spin-torque induced excitations have been demonstrated on uniform 
magnetisations (T), as well magnetic solitons, such as domain walls El, magnetic vortices EH2 and more recently magnetic 
skyrmions I®- Interestingly, the generation of a spin-torque can be via the direct flow of a spin-polarised current perpendicular 
to the plane (CPP) of a magnetic tunnel junction (MTJ) or, as demonstrated more recently, by the production of pure spin 
currents in ferromagnet/heavy metal bilayers due to large spin orbit coupling Emni- An exciting new potential functionality for 
spin-torque nano-oscillators, STNOs, is in the field of high frequency detection, where a spin-torque diode induced d.c. voltage 
due to the mixing of the resistance and the current oscillations im is the result of a radio-frequency (rf) current applied to the 
STNO. This new spintronics-based rectification effect has been shown previously for uniform magnetisation modes in magnetic 
tunnel junctions D3HH, with a sensitivity in excess of the conventional Schottky diode mm and is therefore of significant 
interest as a new type of integrated rf detector or a new generation of STNO read heads of Hard Disks. In this report, we 
focus on spin-torque induced resonant excitations of a single magnetic vortex confined in an MTJ. Indeed, the excitation of the 
fundamental mode of the magnetic vortex, i.e. the gyrotropic motion of the vortex core, has raised a lot of interest in the last 
decade E]0[6j|l5]| where the generation of a steady state excitation of this mode has been achieved by the compensation of the 
damping via a d.c. current. In a more conventional approach, the vortex core can also be resonantly excited via an rf magnetic 
field Q If6]4l8l l or a rf current lH9l if the frequency of the induced core displacement resonantly matches that of the gyrotropic 
mode. An interesting peculiarity of the vortex-based STNOs is the ability to manipulate the chirality and polarity of the vortex 
ED, leading to the promising concept of vortex memories GOl as well as the concept of arrays of magnetic vortices used as 

magnonic crystals ED 1221- 

Vortex core excitation via a spin-polarised rf current 

One of the main purposes of this report is to present a comprehensive study of the effect of a time varying current on the 
magnetic vortex in order to demonstrate the exciting possibilities that exist for STNOs as a new type of nanoscale rf detector. 
Towards this end, we present, in Fig. |T| the rectified d.c. voltage which has been measured as a function of the incoming rf 
frequency. The magnetic tunnel junctions (diameter 500nm) we studied have a magnetic vortex as the magnetic ground state 
(at zero in-plane magnetic field) of the NiFe free layer with a CoFe/Ru/CoFeB synthetic antiferromagnet as the polariser layer 
(see Methods for details). This measurement has been performed for two rf currents, I r f = 0.4 mA (-22 dBm) in Fig. [T] a), and 
Irf = 1.6 mA (-10 dBm) in Fig. [T|b). Both these measurements have been performed with zero applied d.c. current and under 
an applied perpendicular field H perp = 200 mT. For the relatively low rf current, I r | = 0.4 mA, a d.c. voltage, centred at f source = 
110 MHz, is generated when the incident rf current matches that of the gyrotropic mode of the vortex. The lineshape associated 
with this peak in voltage can not be fit with a simple Lorentz/anti-Lorentz peaks as discussed for uniform modes ED, believed 
to be due to the large excitation orbit associated with magnetic vortices and the associated large nonlinearities. The generated 
voltage is due to the spin-diode mixing of the rf current with the TMR induced variable resistance related to the vortex gyrations 
E2- The radius of excitation associated with this peak can be estimated from Ohms law (see supplementary discussion), and 
was found to be p osc ~ 55nm (around 20% of the total radius of the disk). This resonant orbital trajectory of the vortex core can 
also be modelled micromagnetically (see Methods), and the transient trajectory of the vortex core is shown in Fig. |T|c). Upon 
the application of an rf current, whose frequency is that of the gyrotropic mode, f SOU rce = 100 MHz, the vortex core is excited in 
the circular gyrotropic orbit, with a stationary radius of excitation of around 68 nm. When the rf current is increased to I rr = 1.6 
mA, the excitation radius of the vortex core is subsequently increased, which can lead to pronounced differences in the rectified 
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Experimental Simulation Model 



FIG. 1. Vortex core resonance vs. core expulsion. The experimentally measured voltage as a function of the source frequency for a) I,f = 0.4 
mA and b) I,f =1.6 mA, with Idc = 0 mA, and H perp = 200 mT (inset shows voltage response at I t i c = 5 mA). The micromagnetic simulated 
transient trajectory of the vortex core when excited at f sour ce = 100 MHz for c) I r f = 0.4 mA and d) Tf = 1.6 mA. The estimated radius of 
excitation, p, for e) Irf = 0.4 mA and f) Irf = 1.6 mA extracted using Eq.[I] 

voltage response, as shown in Fig. [T]b). When the source frequency is close to the frequency of resonance, 100 MHz < f source < 
150 MHz, the voltage becomes zero in an abrupt transition for this measurement performed at Id c = 0 mA. This abrupt transition 
in voltage can be understood purely in terms of the evolution of the orbit radius of the core. When the rf current is sufficiently 
large, the radius of excitation approaches that of the disk, as demonstrated with the micromagnetic simulations of the transient 
trajectory of the vortex core. Fig. |T]d), which show the radius of the excitation of the core will increase until it reaches the radius 
of the junction, at which point the core is expelled. The system subsequently enters a static quasi-uniform magnetisation state 
(discussed in further detail in Fig. [2]). This vortex core expulsion explains the zero voltage state measured in Fig. |T|b), as there is 
no longer a vortex present, and as such no rf current induced rectified voltage at the gyrotropic frequency. When the frequency 
of the incoming rf current is no longer at the resonant frequency of the vortex core, the free layer returns to the vortex state. Note 
that the vortex expulsion is observable in these junctions because it happens before the vortex core reaches the critical velocity 
for the dynamical reversal of the core J26). 

The variation of d.c. voltage generated by vortex oscillators can be described as AV = AR exp Idc + IrfAR osc cos(ojst). The 
first part is the d.c. voltage variation associated with the change in the mean resistance, AR exp , which occurs when the core is 
expelled. The second part is associated the conventional spin-diode effect due to the mixing of the variable resistance, AR osc 
with the rf current, as demonstrated in Fig. |T|a). Thus, when Id c = 0 mA, the voltage generated by the vortex core expulsion will 
be zero, as shown in Fig. |T]b), however, when Id c 7/ 0 mA, e.g. Id c = 5 mA in the inset in Fig. [T]b), a large voltage response will 
be observed, which is proportional to Id c and highly scalable with AR 0SC and thus with the TMR of the junction. 

Analytical modelling of the resonant core excitation 

As well as comparing the experimentally acquired results with micromagnetic simulations, we aim to further understand the 
observed behaviour by analytically describing the excitations of the vortex core. The force acting on the vortex core can be 
described as F = f’ S | on + f si on + /‘fit + -Poe, where F S j on and A’sion i lf5l l23H25l are the forces associated with Slonczewski 
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torques from the current uniformly polarised in-plane and perpendicular to the plane, respectively, and Jjxt|| is the force due to 
the field-like torque associated with a current polarised in-plane l27l and f'be is the force due to the Oersted field ll28l (further 
details on the individual forces can be found in the supplementary materials). 

For a d.c. current, the force, f ’sionj., acts either with or contrary to the damping force felt by the core, and is responsible for the 
steady state oscillations mentioned previously 130. The forces due to an in-plane polarised d.c. current, F slon j| and F FL t||, will 
result in a static displacement of the core from the centre of the disk, analogous to that of an effective Zeeman field 03, in the 
x and y directions, for Fs\ oa \\ and F FF t respectively (for a d.c. current polarised along The in-plane forces can be expressed 
as fusion| = Asion Ju x and F FL t|| = A|-| i|'./Uy, where Asion and A F lt are the torque efficiencies for the Slonczewski and field-like 


: 0, the ratio of the efficiencies of the two forces can be expressed as 


Asion 

Aflt 


3 b 

2-R£flt 


, where b is 


torques, respectively. At H perp : 

the vortex core radius, R is the disk radius and £ FLT is the relative efficiency of the field-like torque, and has been reported as £ FLT 
= 0.4 for similar asymmetric magnetic tunnel junctions El. The ratio of b/R represents a factor of amplification of the field-like 
torque which is specific to magnetic vortices, due to the field-like torque acting over the body of the vortex, as opposed to the 
Slonczewski torque which will only act only upon the vortex core. The core radius is estimated from micromagnetic simulations 
at b ~ 15 nm, i.e. around twice the exchange length, which results in the ratio b/R = 0.06. As such, the field-like torque will be 
the dominant torque in magnetic vortices compared to the in-plane Slonczewski torque. 

For a time varying current, I(t) = Idc+Irf cos(wsf) (where I r f = J r [7r/i >2 ), the in-plane forces, /* slon and F flx n, will resonantly 
excite the vortex core if the incoming frequency is equivalent to the gyrotropic mode of the vortex core. The radius of these ex¬ 
citations can be found by introducing the rf currents into the Thiele equation (see supplementary materials) and can be expressed 
as: 


Pose 


•Af\/AJ lon + 
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where G is the gyrovector. Mg is the saturation magnetisation, p, is the perpendicular component of the spin polarisation (which 
is proportional to the out-of-plane magnetic field), L is the thickness of the free layer and D is the damping coefficient, J rf and J dc 
are the current densities for respectively the rf and the d.c. current. The total confinement is k(./ ( i c ) = K ms + KocAic, which is the 
sum of the confinement due to the magnetostatic, K ms , and Oersted, no e , terms, u is the excitation frequency and cog = K ' g ' is 
the natural frequency of the vortex. 

Using Eq. |T| we can thus estimate the orbit radius of the vortex core oscillations. For a low current, I r f = 0.4 mA, the radius 
is presented in Fig. |T|e) as a function of the excitation frequency. The maximum excitation radius observed was 63 nm, in good 
agreement with that extracted experimentally (p osc ~ 55 nm). For a larger rf current, our simple model predicts an orbit radius 
larger than the radius of the MTJ (about 260 nm for 1^ = 1.6 mA). This new resonant behavior corresponds to the expulsion of 
the vortex core from the tunnel junction (see in Fig. [T|b)) and can be well described by the analytical model. 

Final magnetization state after core expulsion 

In order to investigate the magnetic configuration of the free layer after the core expulsion, we have studied the influence of 
a small in-plane field, Hip, (obtained by tilting the sample with respect to the perpendicularly applied magnetic field, H pe rp). 
As presented in Fig. [2] a), for H perp = 160 mT, when the angle at which the field is applied is On < 89.6° (0 H > 90.4°), the 
in-plane component of the field is sufficiently strong (Hip ~ 3mT) that the free layer is in a uniform magnetisation state, parallel 
(antiparallel) relative to the polariser. For the angle range 89.6° < 0\\ < 90.4°, the free layer is in the vortex state. 

In Fig. 0 a d.c. current, I dc = 5 mA, is applied simultaneously along with the injected rf current, I rt = 0.8 mA, and the change 
in resistance when the vortex core is expelled is measured. When Hn> < 0, i.e. (t H = 89.65°, and the source frequency is equal 
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to the gyrotropic frequency of the core, the core is expelled to the parallel (P) uniform magnetisation state, i.e. blue dots. The 
opposite is true when Hip > 0, i.e. 0\\ = 90.35°, and the vortex is expelled to the antiparallel (AP) state , i.e. red dots, as shown in 
Fig. |2]b). The ability to tune the final state with a small magnetic field presents the interesting concept of a three state memory, 
where the writing is achieved via a combination of d.c. and rf currents. 

The experimentally observed behaviour is very well reproduced by micromagnetic simulations (I,j c = 4.6 mA and I rf = 0.5 
mA), shown in Fig. |2]c). The simulated micromagnetic final state at three different points are presented, and show that the AP 
and P states (1. and 3.) are in fact C-states where the majority of the magnetisation is directed along the parallel or anti-parallel 
direction. This is consistent with the fact that the resistance does not quite reach the P/AP value represented by the blue/red 
dotted line in Fig. |2]b). 
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FIG. 2. Core expulsion to quasi-uniform state, a) Experimental angular dependence of the resistance showing parallel, vortex and antiparallel 
resistance configurations measured at Id c = 5 mA and H pC rp = 160 mT. b) Experimentally measured frequency dependence of the resistance for 
a range of rf source frequencies. The blue (red) corresponds to the parallel (antiparallel) state, and the green to the vortex state, and the vortex 
state is suppressed when the source is equal to the gyrotropic frequency, f SO urce = 120 MHz. c) The simulated m y component as a function of 
source frequency for an in-plane magnetic field of ±3 mT, where the final m y component is shown for the points labelled 1), 2) and 3). 

The frequency range over which expulsion of the vortex core occurs is asymmetric for the two directions of in-plane field, 
which has been observed experimentally and with the micromagnetic simulations, and can be understood in terms of the shift 
of the vortex core by the in-plane magnetic field being offset by the shift of the vortex core due to the in-plane polarised d.c. 
current induced field-like torque, Fj-ltii (discussed in more detail in the supplementary materials). By varying the d.c. current, 
the frequency bandwidth over which the expulsion is observed can effectively be tuned, and may be potentially interesting in 
terms of accuracy of the rf detector. 

Direct current induced compensation of damping to reduce I™“ 

A crucial feature for potential future applications, is the minimum rf current necessary to observe core expulsion, I™ n , which 
has been measured as a function of the d.c. current at two different magnetic fields, H perp = 0 mT (blue dots in Fig [3]l and H perp 
= 160 mT (red dots in Fig[3j, and the data has been fit by a model based on Eq. [T](see Methods). 

When H perp = 0 mT, the magnetisation of the polarising layer is totally in-plane and as such only the effects of the in-plane 
component of the spin-torque terms need to be considered, i.e. F S ] on || and F FI t . The dependence of I™ 11 is relatively constant 




















6 


as a function of the applied d.c. current. A small negative gradient in the d.c. current dependence of I™ n is predicted by our 
model due to the combination of the d.c. current induced Oersted field altering the confinement felt by the vortex and the force 
due to the in-plane polarised current displacing the core from the centre of the junction (see blue line). 

In the presence of an out-of-plane magnetic field, H perp = 160 mT, the magnetisation of the polarising layer becomes tilted 
out-of-plane, meaning that the out-of-plane Slonczewski term must also be considered, i.e. f ’ S | on | . We emphasise that for the 
Hperp which we applied, the associated spin-torque is not large enough to lead to sustained oscillations of the vortex core. The 
compensation of the damping due to the d.c. current induced Fsion results in a larger resonant excitation radius for a given 
rf current, and as such leads to a strong decrease in I™ n as a function of the d.c. current. This is particularly important for 
applications as it means that the power range over which a working device based on vortex core expulsion can operate can be 
drastically increased with the application of a perpendicularly polarised d.c. current. 



< mA > 


FIG. 3. Reducing I“ m by compensation of damping. The d.c. current dependence of I™" for two magnetic fields, H pelp = 0 and 160 mT and 
0 H = 90.1°, where the data points are the experimental data and the line is the fit with Eq. [T] 

Towards applications: resonance of the core vs. core expulsion 

From the point of view of potential applications as a tunable nanoscale rf detector, an important feature of the core expulsion 
is that the when the core is expelled, the generated voltage will be constant for all values of the rf current, in contrast to 
the relatively smaller orbit spin-diode induced voltage, where the generated voltage is strongly reduced when the incoming rf 
current is reduced. To emphasise this difference, the generated voltage as a function of hr, for a range of d.c. currents, is 
presented in Fig. [4] for H perp = 160 mT. When Ij c = 6 mA, the core is expelled for Irf > 0.15 mA, and the generated voltage is 
constant for all measured values of I rf . The value of I,j c = 6 mA, was the maximum d.c. current that could be be applied before 
degradation of the tunnel barrier. For lower values of Id c , the vortex core is expelled if Irf > I™ n , and displays the characteristic 
constant generated voltage, but where the voltage is reduced due to the dependence on /d c . For smaller rf currents, Irf < I™ n , 
the spin-diode induced voltage is characterised with a sharp decrease as Irf decreases, as the generated voltage will depend upon 
the excitation radius which in turn will depend upon I r f. The current state of the art in terms of generated rectification voltage 
via the spin-torque diode effect is summarised in table [T] The maximum (and minimum) generated voltage, V max (and V rnm ), is 
listed with the corresponding excitation power, P max (and P m j n ). The first two examples are taken from the literature and are for 
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FIG. 4. Radio frequency current dependence. The rectified voltage for observed for the a) core expulsion and b) the spin-diode measurements 
taken at 6 mA and 0 mA, respectively, for a I r f = 0.2 mA. c) The generated voltage, AV, as a function of the excitation current, I,f, for a range of 
d.c. currents, at H pcr p = 150 mT and On = 90.1°. The core expulsion shows a constant AV as a function of Irf, whereas the resonant excitation 
decreases linearly with Irf. 


the in-plane precessional mode HD and the out-of-plane precessional El mode of a uniformly magnetised free layer. The last 
two rows, labelled Resonance (@ 0mA) and Expulsion (@ 6mA), are from this work and represent the best results observed, 
where the values labelled Resonance correspond to the small orbit excitations which decrease with the I r f, as shown in Fig. [4] 
The voltage which is generated from the core expulsion is independent of the excitation power (consistent with Fig. [4]), and can 



Vmax (At V) 

Vrnin (At V) 

P max (dBm) 

P mi „ (dBm) 

Tulapurkar (IPP, MTJ) 07j 

32 

3.9 

-11 

-19 

Miwa (OPP, MTJ) Q3) 

1,300 

2.4 

-26 

-67 

Resonance* (@ 0mA) 

200 

0.63 

-5 

-32 

Expulsion* (@ 6mA) 

2,000 

2,000 

-5 

-43 


TABLE I. Table of current state of the art in voltage rectification measurements available (* = this study) 

clearly be seen to be greater than the previous best achieved in the literature. It is also important to note that this study was 
performed with MTJs which had a relatively low TMR ratio (8.5%), an order of magnitude lower than the junctions presented in 
lU3l (86%). A subsequent further improvement of the maximum generated voltage associated to the vortex core expulsion may 
be anticipated. The maximum effective sensitivity, AV/P r f, possible with the vortex core expulsion in these low TMR MTJs is 
already as large as 40,000 V/W, compared to 12,000 VAV for the OPP modes 03, with another order of magnitude possible 
with tunnel junctions with an enhanced TMR. 

Conclusion 

In summary, we demonstrate the strong interest of the excitation of a magnetic vortex confined in a magnetic tunnel junction 
via a resonant high frequency current both from a fundamental point of view as well as for new approaches for efficient rf 
detection. In particular, we show that, once the generated excitation radius of the vortex core is greater than the radius of the 
junction, vortex core expulsion occurs and the free layer enters the static uniform magnetisation state, leading to a large change 
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in resistance. The rf current required in order to generate an excitation greater than the MTJ radius and therefore observe core 
expulsion can be reduced by displacing the vortex core laterally with an in-plane polarised d.c. current or in-plane magnetic 
field, or by reducing the effective damping by applying a perpendicularly polarised d.c. current. Owing to the highly scalable 
nature of generated voltage associated with the vortex core expulsion, magnetic vortex-based tunnel junctions make a strong 
case as future candidates for a new type of spintronics-based high frequency detector. 

Methods 

The data presented is for circular MTJ nanopillars with a nominal diameter of 500 nm. The stack is composed of // Synthetic 
antiferromagnet (SAF) / MgO (1.075) / NiFe (5) (with thickness in nm) magnetic tunnel junction. The pinned SAF layer is a 
PtMn (15) / CoFe (2.5) / Ru (0.85) / CoFeB (3) multilayer. The NiFe layer is the free layer with a magnetic vortex as the ground 
state. Standard d.c. and high frequency current sources are used to source Id c and I,f, respectively, and the voltage measured 
with a digital voltmeter connected via a bias tee. The perpendicular component of the spin polarisation is estimated at p 2 = 
0.15 at Hperp = 200 mT, from the tunnelling magnetoresistance field dependence. A positive d.c. current polarised in the p z 
direction will act to compensate the damping felt by the vortex core. Before each measurement was made, a large perpendicular 
magnetic field and dc current were applied to reliably initialise the vortex core polarity and chirality, before the frequency of the 
rf source was swept around the gyrotropic frequency of the vortex. The measurements were performed on multiple samples (all 
of similar composition and size) and the same results were obtained Sample 1 (Rp = 28.5 Q) was used for the angle dependence 
presented in Fig. [2] and Sample 2 (Rp = 30.9 fi) was used for the d.c. current and T,| dependence presented in Fig. [4] as well 
as Fig. [] Micromagnetic simulations were carried out with the OOMMF code | [29| . Standard values for NiFe were taken, Mp 
= 580xl0 1 2 3 A/m, Gilbert damping a = 0.009, and the exchange stiffness A = 1.3xl0 -11 J/m. The mesh size was 5x5 nm, for 
a thickness of 5 nm. The polarisation was P = 0.2 and a field-like torque efficiency Git = 0.4. The vortex core was initially 
located in the centre of the 500 nm diameter disk and was excited with a current polarised in the (0 0.8 0.2) direction. The radius 
was modelled by calculating p = p osc + p ( i; sp + p H , and expulsion is assumed when p >0.75R, as observed in micromagnetic 
simulations. The parameters used were Mg = 580xl0 3 A/m, a = 0.009, and A = 1.3x1 0 _11 J/m for a polarisation of 0.2 and a 
Glt = 0.4, which allowed a close correlation between the model and the experiment in terms of generated frequency, and d.c. 
and rf current dependence. 
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